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CN ', Abstract 

o 

£T) . We analyze the multisymplectic Preissman scheme for the KdV 

' equation with the periodic boundary condition and show that the un- 

convergence of the widely-used iterative methods to solve the resulting 
nonlinear algebra system of the Preissman scheme is due to the in- 
_C . troduced potential function. A artificial numerical condition is added 

to the periodic boundary condition. The added boundary condition 
■ makes the numerical implementation of the multisymplectic Preissman 

scheme practical and is proved not to change the numerical solutions of 
the KdV equation. Based on our analysis, we derive some new schemes 
which are not restricted by the artificial boundary condition and more 
efficient than the Preissman scheme because of less computing cost and 
less computer storages. By eliminating the auxiliary variables, we also 
derive two schemes for the KdV equation, one is a 12-point scheme 
and the other is an 8-point scheme. As the byproducts, we present 
two new explicit schemes which are not multisymplectic but still have 
remarkable numerical stable property. Numerical experiments on soli- 
ton collisions are also provided to confirm our conclusion and to show 
the benefits of the multisymplectic schemes with comparison of the 
spectral method and Zabusky-Kruskal scheme. 
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1 Introduction 



The Korteweg-de Veris equation has been used to describe various phenom- 
ena such as acoustic waves in an anharmonic crystal, waves in bubble-liquid 
mixtures, magnetohydrodynamic waves in warm plasma, and ion acoustic 
waves. This equation has two fascinating and significant features. One is 
the existence of permanent wave solutions, including solitary wave solutions, 
and the other is the recurrence of the initial state of the wave form. In 1965 
Zabusky and Kruskal [I] used a finite difference method, i.e. the famous 
Zabusky-Kruskal scheme, to show the existence of solitons which propagate 
with their own velocities, exerting essentially no influence on each other. 
They also discussed the recurrence of an initial state and guessed that the 
KdV equation led to the recurrence. Since then, various methods including 
the finite difference method, the Fourier expansion method 2 and the finite 
element method jH] have been proposed to solve the KdV equation. Unfor- 
tunately, the difference solutions often exhibit nonlinear instabilities when 
a long time integration is performed. In the 1990s, the symplectic schemes 
were introduced and systematically developed for the Hamiltonian systems 
within the framework of symplectic geometry [I]-[n]- Numerical results show 
that symplectic schemes have superior performance, especially in long time 
simulations. The symplectic schemes can be applied to the KdV equation 
which may be transformed into the form of Hamiltonian system. 

Recently, J. E. Marsden etc. and T. J. Bridges etc. [H] proposed 
the concept of multisymplectic PDEs and multisymplectic schemes which 
can be viewed as the generalization of symplectic schemes. Many soliton 
equations such as the KdV equation, the Kadomtsev-Petviashvili equation, 
the Zabolotskaya-Khokhlov equation and the sine-Gordon equation can be 
reformulated into the multisymplectic PDEs and can be solved numerically 
by the multisymplectic schemes. The simplest and basic multisymplectic 
scheme is the Preissman scheme, which has been hot in the last two years 
[9*]-|12|. However, sometimes the direct numerical implementation of the 
Preissman scheme for the multisymplectic equation has puzzled researchers 
all along. When it is applied to solve the periodic boundary problem of the 
soliton equations with degenerate lagrangian like the KdV, K-P equation and 
the water waves equation, the general widely-used iterative method refereed 
as the simple iterative method j2j is not convergent, so are the other iterative 
method such as Newton method and conjugate gradient method. Why? 

To settle the problem, S. Reich eliminated the auxiliary variables to 
get a equivalent scheme for the KdV equation, but he did not consider the 
influence of the boundary condition. 
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In the present paper, taking the KdV equation as an example, we ana- 
lyze the practical computation of the Preissman scheme and find that the 
unconvergence is due to the indeterminacy of the potential function. We 
add a condition on the potential function to fix it up and prove that the 
added condition will not change the numerical solutions of the KdV equa- 
tion. This condition can be stated as a restriction on numerical periodic 
boundary condition. Based on our analysis, we present some new multisym- 
plectic implicit schemes that are equivalent to, but more efficient than the 
Preissman scheme because of less computing cost and less computer stor- 
ages . By converting the implicit term in the multisymplectic schemes to 
an explicit one, we obtain two stable, efficient, explicit schemes for the KdV 
equation. Of cause, they are not multisymplectic any more. An elementary 
but useful method to eliminate the auxiliary variables of the multisymplec- 
tic schemes is also presented and two new multisymplectic schemes for the 
KdV equation are derived. One is a 12-point scheme, the other is an 8-point 
scheme. 

The main purpose of this paper is to develop a method to analyze the 
multisymplectic scheme for the Hamiltonian PDEs and to show how to 
choose the proper numerical boundary condition for the multisymplectic 
schemes and how to derive the new schemes for the PDEs. The method 
presented in the paper can be applied to the Preissman scheme for other 
PDEs and to other multisymplectic schemes. Another aim is to compare 
the performance of the multisymplectic schemes with other kind numeri- 
cal method to see if the multisymplectic schemes benefit the finite difference 
approximations of the PDEs. Thus a series of numerical experiments on soli- 
ton collisions are presented. Compared with the Zabusky-Kruskal scheme 
and the spectral method, the multisymplectic schemes are shown to have 
superior stability, excellent ability to preserve the conservation laws and 
remarkable capacity of long time computing. 

This paper is organized as follows. In section 2 we take a brief review 
of multisymplectic structure of the KdV equation and the multisymplectic 
Preissman scheme. The Preissman scheme is analyzed in section 3, where we 
present an artificial numerical boundary condition for the Preissman scheme 
and verify its rationality. In section 4, some new multisymplectic schemes 
for the KdV equation are derived. Section 5 is for numerical experiments 
and we finish the paper with concluding remarks in section 6. 
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2 Multisymplectic structure of the KdV equation 
and the Preissman scheme 



The general form of the KdV equation with the initial value and the periodic 
boundary condition is 



du du _ 9 d 3 u 

— + + S 2 —^ = 0, t>Q, 

at ox Ox 6 



(2.1) 



u(t = 0, x) = uq(x), u(t, x + a) = u(t, x + b), 

where r] and 5 are two real numbers. 

Introducing the potential <f> x = u, momenta v = 5u x and variable w = 
\4>t + Sv x + V (u), V(u) = r]u 3 /6, the KdV equation (|2.1|) can be rewritten 
as the following Hamiltonian PDEs. 

Mz t + Kz x = VzS(z), (2.2) 
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Each of the two skew-symmetric matrices M and K can be identified 
with a closed two forms. 



x (u, v) =< Mu,v >, 



u 2 (u, v) 



< Ku, V >, 

> is the standard Euclidean 



where u, v are any vectors on R 4 and < • 
inner product on R 4 . 

Both forms to 1 i = 1,2 are closed and therefore pre-symplectic on R 4 , and 
on subspaces where they are non-degenerated, they are symplectic forms. In 
other words, (R 2 ,^ 1 ), and (R 4 ,u; 2 ) are two distinct symplectic manifolds. 
Moreover, each two forms is associated with a different direction, lu 1 is 
associated with time and oj 2 is associated with space. In this sense the first 
order PDEs (J2.2|) is called multisymplectic PDEs or Hamiltonian PDEs. 
The KdV equation is completely characterized by the function S(z), and 
the two skew-symmetric operators M and K. They are all defined on a 
finite dimensional space. 
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The multisymplectic Hamiltonian equation Q2.2|) satisfies the important 
multisymplectic conservation law 



d t [dz A Mdz] + d x [dz A Kdz] = 0, (2.3) 
which, for the KdV equation (|2.1|) . is equivalent to 

d t [# A du] + 2g\ A dw + 5cfo A d«] = 0, (2.4) 

where A is the standard exterior product operator of the differential forms. 

The conservation law Q2.4JI is a strictly local conservation concept that 
does not depend on a specific boundary condition. That is to say, in the 
arbitrary domain of the space-time plane, changes in the wedge product 
dcj)f\du in time are exactly compensated for by changes in the wedge product 
—2(d<p A dw + 5dv A du) in space. 

Bridges and Reich [H] showed that the Preissman scheme for (|2.2|) is 
a multisymplectic scheme which preserves the discrete form of (|2.3|) . The 
Preissman scheme for ()2.2I) is 

\ M KX\ ~ + - = V.S(<$), (2-5) 

where r is the time step, h is the space step, x n , n = 1, 2, • • • , JV; i m , m = 
1, 2, • • • is the regular grids of the integral domain, z™ is an approximation 
t n 7 f T + ^ 7 rn _ 1( m , m\ m +2 _ 1 („m+l i „m\ „ m+ 2_l/_m, 

z™ + i + z™ +1 + ), z = (<f),u,v,w) T , and the corresponding discretized 
multisymplectic conservation law is 

c^+x 1 A (in" 1 ^ 1 - d(p m , Adu m , 



2- 



r 

^Yl 2 A dw™^ 2 + Sdv'n+x 2 A du'^ 2 - d4>n ' 2 A dw' n " ' 2 - ' 2 A dtC 1 2 



,m+i m+i m+= m+ 1 



3 Analysis of the Preissman Scheme 

The Preissman scheme (|2.5|) is an implicit scheme that involves solving a 
nonlinear equations for z m+1 at each time step. The widely-used iterative 
method of this nonlinear equations is as follows. (The analysis and the 
comparison between this iterative technique with other iterative methods 
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such as Newton's Method can be found in El.) 



~{uj +1 + ul+l) -r{w; 



i 



i+i; 2 



+ - 5r(^ - + 2rF'(^), (3.1) 



where i = 1, 2, • • • , n, u{ = \(u{ +1 + u{+l + u{ + u 3 i+1 ), r = ^ is the ratio 
between temporal and spatial steps. 

Here we discuss the numerical boundary conditions. 

By 4> x = u, 



4>(b,t) = <p(a,t) + I u(x,t)dx. 
Set c = u(x, t)dx, then 



— c 
eft 



,6 
J a 

r b 

(cuu x + 5 2 u xxx )dx 



;|u 2 + 5 2 u xx ) x dx 



0. 



We obtain c = f a u${x)dx = a constant. 

Thus the periodic numerical boundary conditions are 
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u j = (u{, u 3 2 , ■ ■ ■ , uiy, v j = (u{,v J 2 ,--- , viy, w 3 = (u{, W2, ■ ■ ■ , wiy, & 

14, • • • , ttf, V(u>\ u^ 1 ) = (Vi, V 2 ,.~, V n f, Vi = V'(l(ui + u> +1 + 
4+1 + ))> ^ = 1,2, • • • V n = V\\{ui + ut 1 + u{ + u{ +1 )), c 

(0, • • • ,0, 2c) T an n dimensional vector, then ()3.1j) and (|3.2I) can be rewritten 
as the following vector form. 

-Au j+1 - B& +1 = --Au j + B& + c, 
2 2 ' 

5Bu j+1 - -Av j+1 = -5Bu j + -Av j , 
2 2 

- 6rBv j+1 + -Aw J+1 - -A<Z> j+1 = 5rBv j - (3.3) 
2 2 

^Aw j - ~A& + 2rV(u J 'X +1 ), 

-Au J+1 + r5w J+1 = -Au j - rBw j . 
2 2 

Also let X = (u J+1 ,v- ?+1 , w- ?+1 ,0 J+1 ) T , we get the iterative form of the 
nonlinear equations 

DX( m ' = 6(X«), Z = 0, 1, - - ■ . (3.4) 

with the coefficient matrix D and the right term 6(X^): (Z), b) = 

( \A -B -\Au 3 +B& + c \ 

-SB \A £Bu-7 - |Av^ 

-5rB \A -\A 5rBv j - %Aw j - \A& + 2tV(u*,uW) 

V \A rB 0, - rBw j ' J 

The initial guess is generally chosen as the value of the previous time step, 
i.e. 



X< > = S«i J .v f .w'.(yj'. 

Unfortunately, this iteration is not convergent because the coefficient 
matrix is degenerated. Note that B is singular, and its rank is n — 1. A is a 
nonsingular matrix only if n is an odd number. To take the row operation 
of the coefficient matrix D, we suppose n be an odd number and A -1 exist. 
Then there is a permutation matrix P, so that 

PD = R, 

where P is the unit lower triangular matrix, R is an upper triangular matrix. 



7 





h a 
2 A 













h a 
2 A 













tA 
2 


I 












The left multiplying (|3,4j) by the permutation matrix P yields an equiv- 
alent system 

DX.V+V = b(X®), Z = 0, 1, • - - , (3.5) 
where the augmented coefficient matrix now becomes (D, b) = P ■ (D, b) = 

-B -%Av$ + B& + c 

- j-BA~ x B ^BA-\B<&> + c) - § Av? 

-| - r(f BA- l ) 2 B r(^BA- l ) 2 (B& + c) - JAw^' - + 2rV 

(3-6) 

Proposition 3.1. The coefficient matrix D of the iteration Q3.4|) is rank 
1 deficient and the deficiency is relative to the potential function <E>. 
Proof By (l3~nT) 

rank(Z)) =rank(D) 

=3rank(,4) + rank([/ nxn + 4<5 2 -^( J B^- 1 ) 3 ] J B) 

Note that A is a full rank matrix and B is a rank 1 deficient matrix, we 
may take a proper r ratio, so that 

det(/ nxn + 45 2 ^(BA- 1 ) 3 ) ^0, 

rank([/ nxn + A5 2 ^(B A'^B) = rank(S). 

Thus D is a rank 1 deficient matrix and it is obvious by the ()3.5|) that 
the deficiency is relative to the potential function <£. 

This proposition implies that we need one and only one more condition 
to fix up the system (|3,4j) and the condition is related to the variable <£. For 
example, 

<p^ + = a constant some a i G [In], (3-7) 

If we add this condition to the (|3.3|) , the rank of the coefficient matrix D is 
just full and the iteration is convergent. It will be clarified by the following 
discussion and the numerical experiments in section 5. 

What we concern about is the solution of the KdV equation (|2.1|) . i.e. 
the variable u J+1 in system ()3.3|) . We hope that the added condition ()3.7|) 
will not change the value of the variable u J+1 in system (|3.3JI . which will be 
proved to be true in the following proposition. 

Proposition 3.2. The variables u J+1 , v J+1 , and B& +1 in system (|3,3|) 
are independent on the added condition (|3.7|) . i.e., the value of <ft{ +1 - They 



S 



are only dependent on , u J , v J . 

Proof. The four iterative equations in the system ((3,5(1 may be written as 
^ u (m) _ = _h AuJ + mj + c (3 8) 

^ v (m) _ ^ 5 ^-i 5 $(i+D = _£ Av J + ^BA-\B& + c), (3.9) 
2 /x 2 /x 

^(m) _ ( ^ + r^A" 1 )^)*^ = r(^SA- 1 ) 2 ( J B^' + C ) 
2 2 n n 

- ^Aw j - + 2rV(uW,u 3 '), (3.10) 

(^U + ^(Sil- 1 ) 3 )^*^ = Au> - ^(BA^f(B^ + c) 

- 4r£A -1 V(u J ', u W ) - |. (3.11) 

We now check the process of the iteration. Recall that the values of u J , 
v- 7 , and B& is given and u^ ) = u- 7 . By (3.11), B<&^ is determined. It is 
independent upon the value of $ + . The determined i?^ 1 ) and (3.8) fix up 
Au^ 1 ). Because A is invertible, uW is fixed up. For the same reason, v^ 1 ) is 
also determined by (3.9). They are all independent upon the value of <p{ +1 ■ 
Substituting u^ 1 ) into (3.11) yields independence of B& 2 \ As the iteration 
goes on, we get three sequences u^ +1 \ v^ +1 ), and B^ l+1 \ I = 0, 1, • • • , 
which are all independent on the value of 0j + . Thus the convergent point 
also has this property, namely, u J+ , v J+1 , and B& +1 are independent on 
the value of . 

This proposition assures the reliability of adding the condition ((3.7|) to 
the system 1)3.3(1 . Actually, the values of u, v and B& are fixed up by 
the Preissman scheme, they have nothing to do with the added condition. 
When to compute the value of the variable we need a condition like ((3.71) 
because the rank of the matrix B is n — 1. For convenience, we may take the 
condition 1)3.7(1 as a boundary condition in practical computation. That is 
<Pi +1 = (or 4^-1 = 0)- The periodic numerical boundary condition ()3.2() 
now becomes 



^=^1 + ^ = 0, (3.12) 



which makes the direct numerical implementation of the Preissman scheme 
practical without changing the numerical solution of the KdV equation. The 
corresponding numerical results on soliton collisions will be presented in the 
next section. Remark: 
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• The above analysis is based on the condition that the number of spatial 
grid points n is an odd number. How to deal with the case with an 
even number? In reference [7j, the numerical experiments presented by 
Marsden et. al. also imply the question: why sometimes the numerical 
results supported on odd spatial grid points are quite different from 
that on even grid points ? 

• By Q3.6JI . we know that the variable w J+1 will change if we change 
the value of <^' +1 in the added condition (|3.7|) . This means that the 
variable w = n</>t + $ v x + V (u) in the Preissman scheme is not fixed 
up. Does the variable w have the important physical meaning? If it 
does, how to fix it up? 

4 Some new multisymplectic schemes for the KdV 
equation 

In this section, we present several new multisymplectic schemes for the KdV 
equation. 

Inspired by (|3.6j) . we obtain a new scheme for the KdV equation (2.1) 



p^ +1 = Mi(q* - \c) + Af 2 (p* + c) + M 3 ( q ' + ^ +1 ) 2 , 
h 4 

q^ 1 = -q> +£(pJ+pJ+l +c ). 
h 



(4.1) 



where p = B&, q = Au, M\ Mi, M3 are three constant matrixs, M\ = 
"2j n + ^r(BA-i) 3 ]~ 1 , M 2 = -^(BA-y - I n , M 3 = -2^M 1 BA" 1 , 

q 2 = {ql,ql--- ,ql) T - 

This scheme is equivalent to the multisymplectic Preissman scheme, so 
it is also a multisymplectic scheme and has an excellent stability. Actually, 
it is composed of the first and the forth line in (|3.6|) . This scheme is more 
efficient than the Preissman scheme because we need not to compute the 
variable w, v and <p. we may take p = and q = Au as new variables, 
furthermore the condition Q3.7JI does not need. It enhances the conclusion 
that u J+1 is independent on the additional numerical boundary value of <fi. 
Moreover, Scheme ( I4.1|) implies a natural iterative form 

p( /+1 ) = Mm? + M 2 (p J + c) + Af 3 ( qJ ^ q( ° ) 2 - Mi^, 

q(W) = _ qJ+ 2 (pJ+p(m)+c)> 
h 
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The computations of the iteration only involves multiplication of matri- 
ces and vectors, It avoids from solving the algebra equations which is the 
main part of computation in other general implicit scheme such as the 
Preissman scheme. After the convergent point (p fe+1 ,q fc+1 ) is obtained, 
solving the equations Au k+1 = q fe+1 yields the numerical solution at the 
k + 1th time step of the KdV equation (|2. 1|) . If we want to solve the system 



B<5> 



k+l 



to get the numerical results of the potential <£ fc+1 , the addi 



tional condition of (p like (|3.7j) is also needed, for the coefficient matrix B is 
rank 1 deficient. 

Eliminating the variable p in the scheme ( 14.1(1 . we have 



h 



,3+1 



+ q J 



= (Mi + ^M 2 )(q'' + q^" 1 ) 



+ M 3 



: s^ )2 + ( i+^ )2 



+ (I + M 2 - -Mi)c. 

(4.2) 



Set = qJ +1 + q», then 



^ +1 = (M 1 + ^M 2 )z^+M 3 



z 3+i 



+ (/ + M 2 --M 1 )c. (4.3) 



Here the matrixs Mi, M 2 and M3 are defined in (|4.1|) , 

This scheme is equivalent to the multisymplectic scheme (|4.1|) , It con- 
tains only the variable u, thus it can be viewed as a multisymplectic scheme 
for the original KdV equation (|2.1j) . Its stability and capacity of long-time 
simulation are the same with the multisymplectic Preissman scheme and 
scheme (|4.1|) . Furthermore this scheme has more benefits such as simple 
form, to practice easily and less computations. 

It is worth mention that even if we modify the scheme (|4.1j) into a real 
explicit scheme which don't need iterations when being applied, the resulting 
scheme still have very nice numerical performance which will be shown in 
next section. The explicit scheme is 



p j+1 = M iq J ' + M 2 p j + M 3 A 2 + (M 2 - \m x )c, 

I h 

qi+l = _tf- + l(pJ + pJ+l + c). 



(4.4) 



Eliminating the variable p yields an explicit scheme for the original KdV 
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equation H2.1|) 



2 v q^ +1 + q*) =(M 1 + -M 2 )(q^' + q^ 1 ) 



+ M 3 



(4.5) 



+ (/ + M 2 --M 1 )c. 



If we modify the implicit term in scheme (|4.2|) into an explicit one, we 
obtain another explicit scheme for the KdV equation 



^(q J+1 + q J ) =(Mi + ^M 2 )(q^' + 



+ M 3 



/tfyj , r q i + q i \ 2 
2 4 J 



(4.6) 



+ (/ + M 2 --M 1 )c. 



But numerical results show that this is an unstable scheme. 

All the schemes above are invalid provided the number n of the spatial 
grid points is even. Next we introduce another method to eliminate the 
auxiliary variables of the multisymplectic Preissman scheme to get two mul- 
tisymplectic schemes for the KdV equation. Both schemes are valid whether 
n is odd or even. 

Let us state the multisymplectic Preissman scheme for equation (2.2) in 
the form 



1 



K-i + K- u i-i- u l ) + -a - K + w i~ w i-i 



i 



2At 
5 



Ax 



1(7 



i-l/ 



0. 



/ j — X I j 7 — 1 



i 1—1 j 



(4.7) 
(4.8) 



+ ^' + ufJ + ~ -^^t 1 + ~ rt-l ~ fi-i) = °> ( 4 - 9 ) 



Ax 



2V' 



.1 



(v 



i-i 

i-X 



j'-i 



(4.10) 
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Eliminating the variable w by Q4.7|) and (4.10), we obtain 

|gK + i - 4+1 + - 4r l ) + -Li - «£f] + ^[^i - Ci 1 

- {4-i - 4-1)) + ^M+i + 44- M + 4r x ) + 4-i + tti 1 ] = 

- 2[V{\{^1 + ^ + 4- 1 + 4)) - v(\(4zi + u\_ x + 4: 1 + 4))}- 

(4.11) 

In the same manner, we may eliminate the variable v by combining (|4.11l) 
and (|4.8|) to obtain 

- 4+1) + 3 K " 4' 1 ) + 3 K-i " «Ti) + u i-2 " "Tal 

+ ^U+i - 4+1) + (*? - 4' 1 ) - (4-1 - 44) - (4-2 - 44)] 

+ ^ i(4+i + 4+1 - M + 4' 1 ) + M-i + 4-1) - (4-2 + 

= -2[y'(i(< 1 1 + + uf 1 + «?)) - V( JccJ + + 44 + «J_ a ))], 

(4.12) 

which together with ()4,9|) yields a new 8-points scheme in only variable u, 
by eliminating the variable (f), 

^1(4+1 + H + H-i + «2- 2 ) - (4+1 + s^r 1 + $4-1 + 
+ [(4+i - H + 3^_! - 4-2) + (<i - H' 1 + 3^ri - «£)] 



+ ^-4 V '(\«4 + «ki + 4' 1 + - v'(\«l + 4-i + C2 1 + 4-2))] 

= 0. (4.13) 



Converting the implicit term in above scheme to a explicit one, we obtain 
a new explicit scheme for the KdV equation whose remarkable numerical 
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property will be shown in next section. 

4^[(<l + 3 4 + H-i + 4-2) ~ (4+1 + Snf 1 + 3t4Zi + <2)1 



5 



2 



+ ^3 [K+i " H + 3ntx - u\_ 2 ) + (uj-± - 3ul L + SulZt - «&)] 
+ ^jy'{\^l + < 1 )) ~ V'i^-l + <2))] = 0. (4.14) 

In the appendix of this paper we present another process of eliminating 
the auxiliary variables to derive a 12-points scheme 



~ 4+1 + 3^ +1 - snr 1 + 3v£l - atCJ + t£] - 

<5 2 



+ aa^^+I - 3 ^ +1 + 3 ^ - 4-1 + M+i - H 

+ n^r 1 + «J + «tJ + 4+d) ~ V'{\(<-1 + 4-2 + 4-1 + 4-i))] 
= 0. (4.15) 

Both schemes (|4.13j) and (|4.15|) are derived from the the Preissman scheme 
(|2.5jl . thus they should be equivalent to each other. Actually they can be 
derived from each other. 

Remark: The method introduced above can be applied to the Preiss- 
man scheme, as well as other multisymplectic scheme, for other Hamiltonian 
PEDs to obtain new schemes. For example, we can get a nine-point scheme 
for the sine-Gordon equation, a six-point scheme for the Schrddinger equa- 
tion, a forty-five-point for the Kadomtsev-Petviashvili equation and so on. 
All these schemes, except for the nine-point scheme, which was discussed by 
Marsden et. al. in are new and expected to have excellent numerical 
stability and capacity of long-time simulation. 



5 Concluding Remarks 

We analyze the multisymplectic scheme for the KdV equation and find that 
the unconvergence of the widely-used iteration method to solve the resulting 
nonlinear algebra system is due to the introduced potential function </>. We 



14 



add a artificial numerical boundary condition on the original periodic nu- 
merical boundary condition. It leads to a new numerical boundary condition 
H3.12|) which makes the implementation of the Preissman scheme practical 
without changing the numerical solution of the KdV equation. The nu- 
merical results obtained with the presented numerical boundary condition 
show the correctness of the condition and the merits of the multisymplec- 
tic schemes. This method for analysis can be easily generalized to other 
multisymplectic schemes and other PDEs. 

We obtain two new multisymplectic schemes for the KdV equation, which 
are equivalent to, but more efficient than the Preissman scheme. We also de- 
velop a method to eliminate the auxiliary variables of the Preissman scheme 
and get two equivalent multisymplectic schemes in only u for the KdV equa- 
tion. One is a 12-point scheme and the other is an 8-point scheme. Com- 
pared with the Zabusky-Kruskal scheme and the spectral method, the new 
multisymplectic schemes are used to simulate the solitary waves. Numeri- 
cal results show that the multisymplecticity do bring the finite differential 
schemes some benefits such as stability, capacity for long time computation, 
and ability to preserve the conservational laws. 

At last we like point out the explicit schemes (|4.5[) and Q4.14JI is also an 
excellent schemes for the KdV equation, they can give the most accurate 
waveforms, which catch well up with those in Figure 6. Furthermore, its 
the stability, capacity for long time computation and efficiency are much 
better than that of the Zabusky-Kruskal scheme, as presented in subsection 
5.5. We are currently analyzing theoretically the stability, conservation and 
other proprieties of the explicit schemes. 

Remark: During the preparation of this paper, Prof. R. MacLachan has 
also derived the same results on the eight-points for the KdV equation. We 
thank him for the discussions and many important suggestions. 
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Appendix 



A The detail process of produce the 12-point scheme 
for the KdV equation 

Let us state the multisymplectic Preissman scheme for equation (2.2) in the 
form 

2A* («<-l + u i ~ u i-l ~ u i~^ + ~k^( W i~ l +w i~ w i-l ~ W i-J = °> 

(A.l) 

^K _1 +4- 4-1 - u Li) - + 4-i + 4' 1 + 4) = °> (a-2) 



~(u> 1 + u\ + u\_\ + u\_ x ) - 1 + ^ - - = 0, (A.3) 

1 . , . . , 1 



" ^(^^ + - 4=1 - vU = 2V'{\{v?r{ + ^ + uf 1 + 

(A.4) 

Taking i = i + 1 in (|A.4|) . we obtain 

-(wi^+wi + w j~l + w 3 i+1 ) - _(^' + <^ +1 - c^p 1 - ^-J) 

(A.5) 

yields 



5 



^ (v^ -1 + 4 + 4+1 + 4+i)) - v'(\(ujzl + <i + 4~' + «f 

(A.6) 
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Let % = i — 1 in ()A.6|) . then 

+Wj~ wlZl ~ Wl_ 2 ) - 2 A ^ Ax (<Pi - V^r 1 - ^-2 + ^2) 

=^ (v'{\{v?r_l + ^ + uf 1 + u{)) - v>{\{v?rl + u{_ 2 + vd + t£i 

(A.7) 

lA~6l+lA~7l 



2 



1 



yields 



7 1 1 7 7 — 1 7 1 7—1 1 J 7—1 7 \ 

+ <+i - <_i - <_! + < + W- - W\_ 2 - W\_ 2 ) 



4Ax v i+l 

- 4AtAz N+i - ^+1 - + ^-i 1 + ^ - ^r 1 - + ^-2> 



^ J— 1 1 i i — 1 i i 7" — 1 1 7—1 . j \ 



2Ax 2V i+1 ^ 1+1 

-t (^(i W" 1 + u * + «£f + " ^'(j + ^-2 + + 4-1 

(A.8) 



and, similiarly, taking i = i + 1 in (|A.2|) . we have 

^^-^-^-^'W+^ + ^-O. <A.») 
calculation of lA^t E Zl leads to 



1 



2Ax 

which implies (i = i — 1) 



:< 1 1 +^ +1 -c 1 1 -^_ 1 )=o. (a.io) 



Ax 



— H" 1 + <4 - 2 »J-i - 2 »m + "CI + < 



' 'v\ 1 + v 3 i -vj_ 2 1 -v J i _ 2 ) = 0. (A.ll) 



2Ax 
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- — ki — - y ields 



^3 {uj + l + u\ +l - Sul -3uj 1 + 3ul_l + 3ul_ x - u\_ 



IK 



3-^ 



i-2, 



1 



(v 3 - , t + ui 



— 



■ 

i-1 



7 i-l 



0. 



In the same manner, taking i = i + 1 in (|A,1|) . we have 

f J ' i 7 7 — 1 7 — 1\ i / 7 — 1 i J 7~ 1 J 

and the sum of the above formula and (1) is 



(A.12) 



0. (A.13) 



n^ 1 



?(»■ 



i-1/ 



0. 



(A.14) 



Combining ifOjl . (fATl . ([A~Ujl . we obtain 



3Ai ( ^ +1 
1 



4+-1 + 3u i - 3u l 1 + 3 ^i-l - 3 «i-l + U \-2 ~ U l-l) 



+ - 3«g - 3u^ + 3^_ x + 3|_ x - 



V?'_2 + Vi-2 1 ) 



\-2 



i-2/ 



Ax 3 m i+1 

=^ (n^r 1 + + + o - v'^ii + ui_ 2 + <z\ + uU) 

(A.15) 

This leads, if j is replaced with j + 1, 
1 



8A ^«i+l u i+i 
1 



+ 3«? +i -3«?+3«?+ 1 -3«?_ 1 + 



4-2 



W 



i-2) 



AAtAx 
5 2 



J+ 1 



i+i 



n^ 1 



Ax 3 
1 

Ax 



K+i + w 



3+1 
i+l 



3uf 



3«? + 3t4_ 1 + 3?Tj 



"i-2 



i-2/ 1 



~ (^'("K + + <i + O - V(i(ii. a + + ^ + t£J) 

(A.16) 
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The sum of the above two formulas is 



3Ai 



1 

AAtAx 



J- 1 



i+i i 3-1 i j-i 



j2 



+ 6«ti + 3^ii + H-\ - H-2 - 4-2 

+ V'(\{< 1 + 4 + <J + - + «J- 2 + ttJ + ttl))) ■ 

(A.17) 

Meanwhile, we take j = j + 1 in (|A.3|) 

+ + 4-i + 4-1) ~^M + - vLi - rfS) = o- ( A - 18 ) 
Sfl yields 

: K +1 + - nr 1 - 4zD - j+j-m* 1 - &l - + = o. 



2At v 4 t ~ 1 2 ^ Ax At 

(A.19) 

Taking % = i + 1 in (jA. 19j) , we obtain 

1 (4 +1 + vRl - 4+1 - 4~ l ) - ^hcMl - - <fcl + vt 1 ) = o- 



2At v 1 1+1 1+1 1 ' Ax At 
By (lA~T9l) + (|A~20l) . we have 



(A.20) 



1 



41} + 24 +1 + 41} - 24- 1 - ut} - vT 1 ^ 



2At +l +1 



AxAt 
which is 

1 



2At 



+ 2u^ + «£g - 2<i X - - 



' - rf" 1 - f^S + rfl}) = 0, (A.22) 



AxAt 
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if i is replaced witlri — 1. 

Combining (|A.21|) . I|A.22|I and l|A.17|l . we obtain a new multisymplectic 
twelve points scheme for the KdV equation (2.1) 

= 0. (A.23) 
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